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Abstract 

We present a study on the properties of electromagnetic plasma waves in 
the region upstream of the Martian bow shock, detected by the magnetometer 
and electron reflectometer (MAG / ER) onboard the Mars Global Surveyor 
(MGS) spacecraft during the period known as Science Phasing Orbits (SPO). 
The frequency of these waves, measured in the MGS reference frame (SC), 
is close to the local proton cyclotron frequency. Minimum variance analy- 
sis (MVA) shows that these 'proton cyclotron frequency' waves (PCWs) are 
characterized - in the SC frame - by a left-hand, elliptical polarization and 
propagate almost parallel to the background magnetic field. They also have 
a small degree of compressibility and an amplitude that decreases with the 
increase of the interplanetary magnetic field (IMF) cone angle and radial dis- 
tance from the planet. The latter result supports the idea that the source of 
these waves is Mars. In addition, we find that these waves are not associated 
with the foreshock and their properties (ellipticity, degree of polarization, 
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direction of propagation) do not depend on the IMF cone angle. Empirical 
evidence and theoretical approaches suggest that most of these observations 
correspond to the ion-ion right hand (RH) mode originating from the pick-up 
of ionized exospheric hydrogen. The left-hand (LH) mode might be present 
in cases where the IMF cone angle is high. PCWs occur in 62 % of the time 
during SPOl subphase, whereas occurrence drops to 8% during SP02. Also, 
SPOl PCWs preserve their characteristics for longer time periods and have 
greater degree of polarization and coherence than those in SP02. We discuss 
these results in the context of possible changes in the pick-up conditions from 
SPOl to SP02, or steady, spatial inhomogeneities in the wave distribution. 
The lack of influence from the Solar Wind's convective electric field upon the 
location of PCWs indicates that, as suggested by recent theoretical results, 
there is no clear relation between the spatial distribution of PCWs and that 
of pick-up ions. 

Keywords: Mars, Upstream, Cyclotron, Waves, Pick-up, Exosphere 
1. Introduction 

Our Solar System is embedded in a plasma flow emanating from the Sun, 
known as the Solar Wind. The Solar Wind expands into the Solar System 
reaching supersonic speeds at a few solar radii from the Sun. The obstacles 
in its path can be classified into three types: absorbers (for example, un- 
magnetized asteroids or the Moon) where the Solar Wind interacts directly 
with their surface, obstacles with intrinsic magnetic fields (such as the Earth 
and Jupiter), and objects with no intrinsic magnetic field and whose atmo- 
spheres directly interact with the Solar Wind. In this group we find planets 
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like Venus and Mars, Saturn's satellite Titan, and also active comets. 

The interaction between the atmosphere of a non-magnetic obstacle and 
the Solar Wind can be described essentially as a non-collisional interaction 
between a magnetic plasma wind and a neutral cloud being ionized by solar 
radiation and charge-exchange with the Solar Wind plasma. An atmo- 
spheric obstacle such as Mars's atmosphe r e gen erates an induced magneto- 



sphere 



A curia et al. 



1998 



Bertucci et al. 



2005| . which has its origin in the 



exchange of energy and momentum between the Solar Wind and the plane- 
tary ions. The induced magnetosphere is preceded by a bow shock because 
of the supersonic nature of the Solar Wind. 

However, the interaction can start far beyon d the bow shock bec ause 
particles from the exosphere (mainly hydrogen) 



Chaufray et al 



2008j are 



ionized several planetary radii away from the object. These particles ar e 



Modolo et al. 



2005|, 



ionized mostly by photoionization and charge exchange 
which add a small amount of energy to the ions with respect to their parent 
neutrals. As the latter are approximately at rest with respect to the planet, 
the ions' planetocentric velocities are also considered to be negligible. These 
newborn ions start gyrating around the interplanetary magnetic field (IMF) 
while preserving their parents neutral's parallel velocity. In the planet's frame 
the gyrating ions also drift perpendicular to B and Ec{Ec = —v x B). 

The physics of planetary ion pick-up is exactly the same as the one seen at 
comets and about which there is a vast literature 



1993) : 



Tsurutani et al. 



1989) : 



Tsurutani 



199l|; 





Vlazel 


e and Neubauer 


Garv 


1991 


|, and references 



therein]. These newborn ions represent a non thermal component of the total 
ion distribution function which is unstable to the generation of electromag- 
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netic waves |Wu and Davidsonl . Il972l : IWu and Hartld . Il974| . The instability 
and wave polarization resulting from the pick up process depends on the IMF 
cone angle av,B, which is the an gle between the solar wind velocity (Vsw) anc . 



the IM F at the time of pick-up 



Tsurutani and Smithl . 



1986 



Tsurutani et al 



1987|. If Vsw is parallel to the IMF, the newborn ions will form a beam in the 
solar wind frame and the electr omagnetic io n-ion right-hand (RH) resonant 



instability will be predominant 



Gary 



1993] ■ On the other hand, if Vsw is 



perpendicular to the IMF, the newborn ion distribution function will drive 
the electromagnetic ion-ion left-hand (LH) mode unstable. Both instabilities 
have maximum growth rates at k x B = 0, where k is the propagation wave- 
number. At moderate angle s of av,B the RH instability is still predominant. 



Brinca and Tsurutani 



19891 found that the maximum growth rate of the LH 



instability is 
according to 



larger than that o: 



Convery and Gary 



the RH instabihty for ay n > 7 5°, whereas 



1997 | 



Gary and Madland 



1988) ■ this cutoff 



Garv et al.. 


1989; 


(1, u 1 

Brinca. 


/|— ^ 

1991 



cone angle is a^.B = 90°. 



is approximately zero for moderate ay^B 
In this expression fij is the newborn ion gyrofrequency, {7^°" is the ion drift 
velocity along the magnetic field -B, and lo and k are respectively the wave 
frequency and wave vector in the Solar Wind frame in the case of propaga- 
tion parallel to 5. As the spacecraft has a negligible planetocentric velocity 
compared to that of the Solar Wind Vsw, we obtain the following expression: 



OJ - k ■ vff] 



■ ^] k 



(1) 



where Usc is the frequency of the plasma wave in the SC frame and k = k/\k\. 
Since the frequency of the electromagnetic plasma wave in the newborn 



ion reference frame is Uion = OJ — k ■ v^°p, and the velocity of the planetary 
particle before ionization is negligible with respect to Vgw, we obtain the 
useful expression: 

Use = -^i (2) 

where Qi = {qiB/nii), with B = \B\ , and q and m being the charge and mass 
of the newborn ion respectively. This means that the waves generated by the 
RH instability are characterized by a frequency which, in the SC frame, is 
close to the local ion cyclotron frequency with a left-handed polarization. 

When the parallel propagating ion-ion LH instability prevails, the Doppler 
correction is very small and therefore the waves will be observed with a left- 
hand polarization in both the plasma and SC frames at the local ion cyclotron 
frequency. 

This suggests that the occurrence of waves at the local ion cyclotron 
frequency of a particular ion species in the SC frame can be associated with 
the occurrence of the pick up of such ions. In this sense, the presence of 
plasma waves is, a priori, an important diagnostic tool for the evidence of 
ionized exospheric particles. 



The first observation of 
by Phobos-2 [Russell et al. 



'CW s upstream from Mars' bow shock was made 



1990|. These waves had small amplitudes (~ 0.15 



nT), they were left-hand elliptically polarized in the SC frame, and propa- 
gated at a small angle to the mean magn etic field. P CWs have also been 
observed by Mars Global Surveyor (MGS) 



Brain et al 



2002 



Mazelle et al 



2004| . The frequency, polarization and propagation angle of the waves de- 



tected by MGS were similar to those determined from Phobos-2 observations, 
except their amplitude (2-3 times greater). Waves at the local proton cy- 



5 



clotron frequency have al so been observed upstream from Venus 



2009| and active comets 



Tsurutani 



199l|; 



Delva et al. 



Mazelle and Neubauei 



1993) and 



references therein]. 

The eccentric orbits of MGS during the mission's pre-mapping phase 
2OO1I allowed obs e rvatio ns of PCWs up to 15 Rm (li?A/=3390 



Albee et al. 



km: Mars radius). iBrain et al.l 2002j performed a statistical analysis of the 
properties of these waves during the first aer obraking (ABl) and scie nce phas- 



Wei and Russell 



2006l | analyzed 



ing orbit (SPO) phases. A few years later 
85 events during the ABl phase and discussed the generation mechanisms at 
large distances, as well as the possible distribution of waves depei iding on the 



direction of E^. = —v x B. Similar analyses were presented by 



Delva et al 



2009|, based on measurements provided by the magnetometer onboard the 
Venus Express spacecraft during 450 orbits. 

In this study we carry out an analysis of the PCWs detected by the 
magnetometer and the electron refiectometer (MAG/ER) onboard MGS in 
the region upstream from the Martian bow shock during the mission's SPO 
phase. We analyze the frequency, propagation and polarization properties of 
these waves and discuss the generation mechanisms and their relationship to 
the neutral densities at the exosphere of Mars. We also analyze the spatial 
distribution of these waves in a magneto-electric coordinate system centered 
on Mars (MBE). Finally, we study the im plications of our results a nd compare 



;hem with recent studies around Mars 



Wei and Russell 



Delva et al. 



2006) and Venus 



201l|. The article is structured as follows. In section (2) we 
describe the capabilities and limitations of MAG/ER in characterizing the 
properties of the PCW's, along with a description of the various methods of 
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analysis that we applied to the measurements. In section (3) we show typical 
examples of PCW's as detected by MAG/ER and obtain their properties 
based on the methods described in section (2). Next, we show statistical 
analyses of the amplitude and the spatial distribution of these waves, as well 
as the Solar Wind IMF cone angle associated with them. In section (4) we 
present a discussion of the results, and the theoretical approaches that might 
explain the generation of these waves. Finally in section (5) we summarize 
our conclusions. 



2. Upstream waves: Analysis Methods 



MGS entered into orbit around Mars on 11 September 1997 



Albee et al. 



2001| . During the pre-mapping AB and SPO phases, MGS provided measure- 
ments of the Martian environment from the unperturbed Solar Wind down 
to the neutral atmosphere from 1683 elliptical orbits. After these orbital 
phases, MGS reached a final circular mapping orbit at 400 km altitude. 

MGS carried a combination of a twin-triaxial fiuxgate magnetometer 
system (MAG) and an electron spectrometer used as a refiectometer (ER) 



Acunaetal 



1992| . MGS did not carry any instruments dedicated to the 



measurements of ion properties. The magnetometers (MAG) provided fast 
measurements (32 vectors/s) over a wide dynamic range (from ±4 nT to 
±65536 nT), and the electron spectrometer (ER) measured the electron fiuxes 
in 30 logarithmically spaced energy channels r anging from 10 eV to 20 keV 



Mitchell et al. 



200l|. In this 



with a maximum integration resolution of 2 s 
study we analyze the averaged fluxes over all directions with energies higher 
than 30 eV (below this energy value, the electron distribution function is 
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affected by spacecraft pfiotoelectrons) . 

Tfie complete removal of spacecraft fields from the upstream MAG mea- 
surements was difficult, as the sensors were not mounted on a boom, but 
on the outer edge of each of the two solar panels. However, s pacecraft fielc 
mode ling allowed the reduction of their influence to ±1 nT 



AcunaetaL 



2001 |. 



We used the MAG/ER measurements from 27 March 1998 through 24 
September 1998. During this period, MGS performed 372 elliptical orbits 
grouped into two subphases: 

• From P202 (27/03/98) to P327 (27/05/98). [SPOl sub-phase] 

• From P328 (27/05/98) to P573 (24/09/98). [SP02 sub-phase] 

The P prefix indicates the periapsis number. 

The SPO orbits had a constant period of 11.6 hours with apoapses above 
the south pole at distances of roughly 6 Rm- They also displayed a mono- 
tonic change in their local time from noon to 4AM. During the same period, 
the Martian season varied from early southern hemisphere summer to mid 
northern hemisphere spring. The SPOl and SP02 sub-phases arise as a re- 
sult of the presence of a hiatus due to a solar conjunction around July, 1998. 
As a result, SPOl orbits had local times between noon and 10 AM, and 
SP02 orbits had local times between 10 AM and 4AM. Also, SPOl orbits 
are characterized by low to moderate solar zenith angles (SZA), whereas up- 
stream observations during SP02 were closer to the terminator plane (see 



figure l,lBrainelalJ |2Q02|). 
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The sampling frequency of MAG measurements analyzed here is fsamp = 
1.33Hz, and due to the influence of spacecraft fields, the study is focused on 
periods where the strength of the IMF is equal to or larger than 4 nT. 

Analysis 

The analyses performed on MAG/ER data consisted in a characterization 
of their spectral properties, their polarization, their degree of coherence, their 
magnetic connection to Mars's bow shock, and their spatial distribution in 
an electromagnetic coordinate system. A more detailed description of them 
is given below. 

2.1. Dynamic spectra and correlation 

We generate dynamic Fourier spectra of the magnetic field components in 
the Mars Solar Orbital (MSO) coordinate system. This coordinate system is 
centered at Mars with its Xmso axis pointing toward the Sun, Zmso being 
perpendicular to Mars's orbital plane and positive to the ecliptic north, and 
while Ymso completing the right-hand system. 

We also calculate the cross correlation between the fluctuations present 
in the electron flux measurements for specific energy channels, and in the 
component of the magnetic field along its mean value. 

2.2. Polarization - MVA 

The wave- vector and polarization of PCW's were obtained from minimum 



variance analysis (MVA) 



Smith and Tsurutanil . 



1976 



Khrabrov and Sonnerup 



1998| . This method provides an estimate of the direction of propagation for 



a plane wave by calculating the eigenvalues of the covariance matrix of the 
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magnetic field within each interval ( the maximum, intermediate and mini- 
mum eigenvalues are Ai, A2 and A3, respectively). Then, the wave- vector k is 
associated with the minimum variance eigenvector. Additionally, the angle 
between k and the mean magnetic field were calculate d with an error based on 
the n umber of measurements and eigenvalues ratios 
1998 |. 



Sonnerup and Scheibld . 



2.3. Polarization and coherence - Coherence Matrix 

An additional tool of analysis consisted in calculating the coherence ma- 
trix of the magnetic field for a 0.015 Hz interval around the local proton cy- 



clotro n frequency (MAG uncertainty) 



Fowler et al 



19671 : 



McPherron et al 



972|. The elements o: 



Tsurutani et al. 



this matrix provide information about the coherence 
2009|, polarization and ellipticity within the selected fre- 
quency range. These parar neters are explicitly indic ated in equations (21), 
(20) and (23), respectively [Rankin and Kurtzl . Il970 |. 



2.4. Magnetic connectivity with the Martian bow shock 

We also determine the spatial location of wave observations with respect 
to the Martian foreshock. For this purpose, we use a static model of the bow 



shock 



Vignes et al. 



2000( and look for an intersection point with the field 



line sampled by the spacecraft, assuming a uniform IMF. 

2.5. Spatial distribution 

Several theoretical studies 



Modolo et al 



2005| suggest a dependence of 



the spatial distribution of pick-up ions upon convective electric field E^. We 
investigate if a similar pattern is found on the distribution of waves. As- 
suming Vsw = —400 km/s Xmso, we study the spatial distribution of waves 
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(including their amplitude) , with respect to Ec, by introducing a "electro- 
magnetic" coordinate system (MBE) which is centered at Mars and where 
the Zmbe axis is parallel to Ec-, Xmbe is antiparallel to Vsw, and Ymbe 
completes the right hand triad. The IMF cone angle is also included in this 
analysis. 

These methods have been applied on selected time intervals to illustrate 
the main properties of these waves and then on 10-minute overlapping inter- 
vals in order to study the changes in the waves' properties along the MGS 
trajectory. The duration of these segments was chosen so as to contain a 
significant number of proton cyclotron periods. 

3. Results 

Figures [Hand [2] show typical examples of wave activity present in MAG/ER 
measurements. Figure [1] illustrates an example of PCW's seen by MAG up- 
stream from the Martian bow shock. The top three panels show the magnetic 
field components in the MSO reference frame. At first glance, the oscilla- 
tions show a well-defined frequency, large amplitudes (up to 4 nT) and a 
high degree of coherence. In this interval, the average magnetic field (Bo) is 
quite steady: its magnitude is slightly above 8 nT and its orientation makes 
an angle of 34°with respect to the Mars— Sun direction. Figure [2] shows the 
electron fiuxes measurements seen by ER for energies 116, 191 and 314 eV 
(with maximum time resolution of 2s). Oscillations with features similar to 
those measured by MAG can be observed, which we assume also corresponds 
to fiuctuations in the total electron density. 
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Figure 1: Magnetic field measurements during part of tire orbit P216 in MSO coordinates, 
(April 3, 1998). 



3.1. Temporal spectra 

When simultaneous MAG and ER measurements are compared, we find 
that the component of the magnetic field parallel to B^, and the omnidi- 
rectional electron flux are in phase. Figure [3] shows the fluctuations in the 
parallel component of the field (upper panel), and in the electron flux at 116 
eV (lower panel) during part of the orbit P232. The cross-correlation of the 
two time series shown in the enclosed panel displays a peak at zero displace- 
ment with a correlation value of 0.7, falling down to 0.44 at ±2s. We repeat 
this analysis for other events and for different energy channels with similar 
results, which is consistent with both time series being in phase, considering 
an error of ±2s associated with the cadence of the ER instrument. 

We also find that most of the waves observed by MGS MAG/ ER have 
frequencies in the SC frame which coincide with the local proton cyclotron 
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Figure 2: Electron flux measurements during part of the orbit P233, (April 11, 1998). 



frequency {Qp = B/mp), within a ±1 nT uncertainty. Figure S] shows a 
dynamic Fourier spectrum of the Ymso component of the magnetic field 
(Bymso) in the orbit P216. The peak in the frequency of the oscillations is 
systematically close to the calculated proton cyclotron frequency for several 
hours. 

3.2. Polarization - MVA 

MVA shows that the waves are planar, propagate almost parallel to the 
background magnetic field and have a left-hand circular polarization (in the 
SC frame) with respect to B^,. Figure [5] shows the results of MVA applied on 
10 cyclotron periods of data during orbit P216. In this case, MVA yields a 
high A2/A3 ratio (50.3) and a small angle between k and B^ {O^b = 8°). The 
hodogram describing the maximum and intermediate variance components 
clearly shows a left hand polarization, as the mean magnetic field for this 
interval points out of the page. 
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Figure 3: Fluctuations in the parallel component of the magnetic field and in the electron 
flux measurements (E=f 16 eV) during part of the orbit P232, (April 11, 1998). 



In the case of orbit P216, the propagation properties of the waves remain 
unchanged for several hours. As a result, an estimatio n of the 6kR from 10 



Knetter et al. 



2004 1 



proton cyclotron periods windows for which A2/A3 > 10 
yields a value of 9.9°. 

We also find that these waves typically have more power perpendicular 
to the mean- field direction than along it, which again indicates that k makes 
a small or moderate angle with the mean-field direction. 
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Figure 4: Fourier dynamic spectrum of Bymso, Orbit P216, (April 3, 1998). The 
calculated local proton cyclotron frequency from MAG data is plotted in black for reference 
(line in the middle). The black lines at the sides correspond to the error bars associated 
with the MAG data. P is the power spectral density. 



3.3. Polarization and coherence - Coherence Matrix 

The apparent coherence of the waves shown in Figure [1] is confirmed with 
an estimate of the coherence coefficient described in section (2). Figure [6] 
illustrates an example of PCW's (orbit P204) with a high degree of coher- 
ence (> 0.7) and polarization (> 70%) for over an hour. Their ellipticity 
coefficient indicates a left-hand elliptical polarization (< -0.6). Similar anal- 
yses applied to other SPOl orbits suggest that these properties also remain 
unchanged for long periods of time (typically 1 hour). 

We also see that the wave amplitudes decrease with radial distance from 
the planet, supporting the idea that Mars is the source of these waves. Figure 
[7] illustrates what is observed in orbit P216. A similar behavior is observed 
for other orbits. 
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Figure 5: MVA hodogram corresponding to a 10 cyclotron period interval during orbit 
P216. The mean magnetic field points out of the page {Bo = [-1.02,-0.18,8.27] nT). 
The cross indicates the start of the time series analyzed. 

3.4- IMF dependence: Connection to the Martian bow shock 

The properties of PCW's are the same inside and outside the nominal 
Martian foreshock. The results are summarized in Table 1 which shows the 
average properties of PCWs for selected orbits that were chosen based on 
their particularly high values of degree of polarization and coherence. For 
example, orbits P204 and P207 have similar wave properties: they both 
are characterized by a high degree of polarization (> 90%), a large coherence 
(0.9) and a left-hand ellipticity (-0.72). 6kB is less than 15°and their IMF cone 
angles differ in approximately 8°. Furthermore, in both cases the waves are 
observed in the same region (approximately between (2.09,-1.27,-2.22) Rm 
and (1.7,-1.38,-4.21) Rm in MSO coordinates). However, whereas in P204 
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Figure 6: Properties of the ultra low frequency waves observed by the MGS spacecraft 
during a part of the orbit P204, (March 28, 1998). 



the magnetic field lines are connected to the bow shock, this does not occur 
in P207. The difference arises because in the first case. Bo is approximately 
(3,-1.5,-3) nT and in the second one (3,-5,0.7) nT. This strongly suggests 
that these are not associated with Solar Wind backstreaming ions within the 
ion Martian foreshock. 

3.5. Overall occurrence 

A study of the overall occurrence demanded the definition of a few criteria 
of the type of oscillation that would be considered to be a PCW. These 
criteria are: 

• Frequency close to the calculated proton cyclotron frequency (±lnT). 

• MVA: events with A2/A3 > 10. 
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Figure 7: The relative wave amplitudes decrease with radial distance from the planet. 
Orbit P216, (April 3, 1998). An exponential fit y = Aexp{—^) yields a scale height 
a = (4.15 ±1.5) i?M- 



• Events with coherence> 0.7, ellipticity < —0.5 , degree of polarization 
> 70%. 

As we explained in section (2), we study the temporal variation of the 
different waves properties breaking time series into 10- minute overlapping 
segments. 

Figure [8] shows the wave amplitude versus planetocentric distance for all 
SPO orbits. The amplitude is estimated from the power spectral density in 
an interval centered at the proton cyclotron frequency with a width of 0.015 
Hz. Cyan crosses correspond to SPOl, while brown open circles indicate 
wave occurrence for SP02. The black curve displays the average amplitudes 
for 0.5 -Rm bins, showing a general decrease of amplitude with increasing 
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distance. Assuming that the source for these waves are exospheric pick-up 
protons, the obtained result is to be expected for at least two reasons: the 
ion density decreases with distance, and for smaller MGS altitudes, waves 
have more time to grow as they are convected by the Solar Wind. 
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Figure 8: Amplitude of the PCW's as a function of the ahitude for all SPO orbits. 
(Criteria stated at the beginning of section 1575)) . 

Figure [9] shows the amplitude of the waves (estimated in the same way as 
in figure [S]) as a function of the observed IMF cone angle for all SPO orbits. 
The cyan crosses correspond to SPOl and the brown open circles to SP02, 
while the black curve displays the average amplitudes for 10° bins. This figure 
does not show any statistical shift between these two populations. The 10° 
bin curve shows a slight shift toward small amplitudes as the IMF cone angle 
increases. This trend is similar to the behavior of the saturation amplitude 
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of the instablity shown in ICowee et aLl 2012| . However, ICowee et aLl 2012 ] 
suggested that the saturation level is likely not reached in the upstream 
region of Mars, which is not necessarily the case for the waves measured by 
MGS. The observed differences could be due to the non homoscedasticity 
of the av,B distribution and the ay^B estimation errors, which make direct 
comparison with their results not so straightforward. 
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Figure 9: Amplitude of the PCW's as a function of the IMF cone angle for all SPO orbits. 
(Criteria stated at the beginning of section [ 



We also analyze the spatial distribution of the orbit segments where waves 
were found, making use of the MBE coordinate system. Figure [THla shows 
the wave amplitudes (color coded) for the projection of the trajectory on the 
{yMBE, Zmbe) plane for all SPOl and SP02 orbits. Figure [TOlb shows the 
av,B value instead. In the absence of velocity vector measurements, av,B is 
obtained assuming that Vsw is parallel to the Sun-Mars line and therefore: 
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Oiv,B — acos{ ^^^° ). In both figures, closed and open circles correspond to 
SPOl and SP02 respectively. The analyses do not show any clear difference 
in the spatial distribution of the waves between these two sub-phases: 60% of 
the waves occur above the Zmbe — plane and 40% are below it, while the 
spacecraft is 48% and 52% of the time in each hemisphere. As a result, the 
spatial distribution of PCWs does not seem to be affected by the Solar Wind's 
convective electric field. In fact, the spatial occurrence of the amplitude of 
the waves and the cone angle av,B associated with them are not affected by 
the electric field either. Moreover, these waves are observed even when Ec is 
relatively weak. 

Another interesting aspect of this study is the clear difference in the 
occurrence rate of the waves between SPOl and SP02: whereas PCWs are 
present in 62% of the SPOl upstream observations, they only appear in 8% 
of the time spent by MGS in the upstream region during SP02. The waves 
during SPOl last for longer and have a higher degree of polarization and 
coherence than those observed in SP02. 

This difference in the occurrence of PCWs is in principle not due to 
changes in the av,B as the distributions for both subphascs peak at around 
60° (Parker's spiral angle at Mars 55°) with a standard deviation of 19°. 

4. Discussion - Planetary Ion Pick-up and Waves at the Cyclotron 
Frequency 

One of the main results of the present study is that the waves in the 
region upstream of the Martian bow shock are observed, in the SC frame, at 
the local proton cyclotron frequency and are left-hand elliptically polarized. 
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Figure 10: Trajectory of the spacecraft SPO orbits with PCW's (MBE coordinate frame). 
Upper panel (a): Color coded by the amphtude of the waves. Lower panel (b): Coded 
by color by the cone angle associated with the MAG measurements where it is observed 
PCW's. Closed and open circles correspond to SPOl and SP02, respectively. 
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As described in the introduction section, the Vlasov-Maxwell theory pro- 
vides a mechanism capable of generating this type of waves as a result of 
the interaction between the magnetic field and t he distribution function of 
"exospheric " pick-up protons in the Solar Wind 
particular, numerical studies 



Russell et al. 



Gary 



1990|. In 



1993| suggest two possible modes LH 



and RH whose predominance depends on the initial distribution function of 
the pick-up protons, which in turn, depends on the IMF cone angle. The 
dispersion obtained in the values of the IMF cone angle (where Vsw was 
assumed to be along the Sun-Mars direction) indicates that in principle both 
modes could be present in MAG observations. On the other hand, all wave 
observations in the SC frame displays a LH polarization. In the absence of 
ion measurements it is impossible to obtain a value for the phase velocity 
of these waves. In particular, the angle between k and Vsw (that is neces- 
sary to compute the Doppler correction) at which the p olarization changes 



Cowee et al. 



2012 | 



remains unknown. However, following theoretical results 
one could assume that the phase speed is of the order of the 10% of the Vsw, 
in which case the cutoff angle would be around 95°. In that case, most of the 
waves observed would correspond to the RH mode. LH waves have been ob- 
served a t comet Gr i gg-Sk jellerup 



near 90° 



Cao et al. 



Neubauer et al. 



2003j for IMF cone angles 



1998| and references therein]. 



MVA analysis shows that the observed PCW's are planar and propagate 
almost parallel (but not exactly) to B. This could be the result of nonlinear 
effects not addressed in this paper. This non zero 6kB is responsible for a 
small compressibility (typically, 6Bpar/B ~ 0.25). 

The compressibility of these waves was studied from the cross-correlation 
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between the fluctuating electron flux (with the electron fluxes taken as a 
proxy of the electron density) and the parallel magnetic field component. The 
maximum values for the cross-correlation coefficient between both time series 
are associated with phase shifts between and 2s. The quality of this result 
is limited by the following two factors: the influence of the magnetic fields 
generated by the spacecraft and the time resolution of the ER instrument (2s). 
These two sources of error impose a tradeoff in our aim to obtain the best 
possible cross-correlation between the two time series. With this in mind 
we calculated the cross-correlation for several intervals with mean magnetic 
fields between 4 and 8 nT. For all the cases considered, the results indicate a 
maximum in the correlation coefficient located within 2s from zero time lag. 
This means that there is a systematic correlation between the two data sets 
within the errors caused by the spacecraft fields and the ER time resolution. 
But also, these results indicate that anticorrelation is not possible (errors are 
not as high as half a proton cyclotron period). 

The analysis of MAG measurements also shows that the frequency, polar- 
ization and coherence of these waves do not seem to depend on the direction 
of the IMF (see Table 1). It is then not surprising that the observed waves 
are not associated with the foreshock. In addition, the amplitude of the 
observed waves shows a decrease as the IMF cone angle increases. 

The lack of influence of Ec on the spatial distribution of PCWs as well 



as their amplitudes suggests that the link betwe en the spatia 



Cowee et al. 



distribution of 



2012) . On the 



the pick up ions and the waves is not obvious 
other hand, the decrease of the observed wave amplitude with radial distance 
from the planet, supports the idea that Mars is indeed the source of these 
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waves. 

We also analyzed the PCWs overall occurrence using very strict criteria. 
We found a very clear difference in the percentage and the properties of 
PCWs between SPOl and SP02. However, a statistical study of the IMF 
cone angle did not show any significant difference between both periods, 
suggesting that the pick-up geometry might have not changed significantly 
from SPOl to SP02. The difference in the occurrence of the waves between 
SPOl and SP02 might be related to temporal or spatial changes in the 
pr operties of the Martian hydrogen exosphere. 



Wei and Russell 



2006j studied the wave occurrence in ABl, as seen from 



a reference frame that takes into account the convective electric field. They 
observed that the PCWs at large distances (85 events present in 9 orbits 
with B > 5.6 nT) occur intermittently and predominantly in the hemi- 
sphere where Ec points out from the planet (+E hemisphere). In order to 
explain this behavior the authors proposed a mechanism where an exospheric 
hydrogen atom is ionized and then accelerated by the Solar Wind electric 
field. A charge exchange collision would then produce fast neutrals able to 
reach distant regions where they are re-ionized and generate cyclotron waves 
downstream of the planet and on the + E hemisphere. 

During SPOl and SP02 the wave spatial distribution for distances smaller 
than 6 Rm does not seem to depend on the orientation of the Solar Wind 
convective electric field (even after more restrictive criteria have been ap- 
plied). Furthermore, waves are found even when Ec is relatively weak. We 
also observe PCWs far away into the -E hemisphere, although there is no 
known mechanism to move ions against the electric field. One explanation 
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to this is that the wave distribution does not necessarily follow the one of 
pick-up ions. It is important to note that pick-up ions generate waves whose 
wavelengths are of planetary scale. Ion density is not likely to be homoge- 



neous over such large spatia . 



obtained in 



Wei and Russell 



distances. The difference between the results 



2006| and the ones presented here could be due 



to the fact that MGS sampled different regions during ABl and SPO. On 
the other hand it is also important to point out that the lack of correlation 



between the spatial distribution o: 
has also been observed at Venus 



the PCWs and the convective electric field 



Delvaetal 



2011 | 



5. Conclusions 

In this study we analyzed upstream MGS MAG/ER measurements during 
the pre-mapping SPO phase. Analyses of these measurements show waves 
whose frequencies in the reference system of the satellite are near the local 
cyclotron frequency of protons. These ultra low frequency plasma waves are 
also characterized (in the spacecraft frame) by a left-hand elliptical polar- 
ization and by an amplitude that decreases with the increase of the IMF 
cone angle and the radial distance from the planet. They propagate almost 
parallel to the background magnetic field and show a small degree of com- 
pressibility. The observed waves are not associated to the foreshock and their 
properties do not depend on the angle a^B between the solar wind velocity 
and the magnetic field direction. 

The most plausible explanation for the existence of these waves is that 
they are generated from the pick-up of newborn protons from the Martian 
exosphere whose distribution function kinetically interacts with the Solar 
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Wind plasma and the interplanetary magnetic field [Russell et al.l . Il990 |. 

From a theoretical perspective, there are two possible modes that can be 
generated in the plasma frame which depends on the pick-up proton distri- 
bution function. For relative orientations between Vsw and B ranging from 
parallel to almost perpendicular, the RH mode will predominate, while for 
almost perpendicular orientations the LH mode will dominate. The reason 
why we always observe left-hand polarized waves in the SC frame is because 
of their corresponding Doppler corrections. 

The spatial distribution of these waves does not seem to depend on the 
orientation of the Solar Wind convective electric field and they are also ob- 
served even when it is relatively weak. PCWs are also found far into the 
negative convectiv e electric field regio n. This supports the idea from recent 



simulation results 



Cowee et al. 



2012| that, even though the origin of these 



waves are the pick-up ions, there is no clear relation between their corre- 
sponding spatial distributions. 

We also found a clear difference in the PCWs occurrence between SPOl 
and SP02. Moreover, it is observed that the waves observed during SPOl 
last longer, and have higher degree of polarization and coherence than those 
observed in SP02. Since there is no significant difference in the IMF cone 
angle distribution between both periods, these differences might be due to 
changes in the properties of the exosphere and/or ionization rates. 

A more accurate determination of the wave modes present at Mars will 
require information about the pick up ion distribution function as well as 
that of the Solar Wind. We hope MAVEN spacecraft (Mars Atmosphere and 
Volatile Evolution Mission) will provide such measurements contributing to a 
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better understanding of the microscopic processes arising from the interaction 
between the Solar Wind and Mars. 
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Orbit 


D. Pol 


Ell 


Coh 


A2/ A3 




|B| 


OiV,B (°) 


Connection 


204 


94.7 


-0.72 


0.90 


31.4 


8.9 


4.86 


51.9 


Yes 


207 


95.4 


-0.72 


0.90 


32.8 


14.2 


5.94 


59.6 


No 


208 


89.6 


-0.67 


0.80 


17.4 


10.0 


5.04 


54.1 


Yes 


209 


88.6 


-0.67 


0.78 


13.9 


19.3 


9.04 


49.8 


Yes 


211 


91.7 


-0.68 


0.83 


15.3 


15.6 


4.65 


56.6 


Both 


215 


85.4 


-0.65 


0.73 


18.1 


38.4 


13.5 


40.2 


No 


216 


94.5 


-0.69 


0.89 


20.2 


17.6 


7.31 


33.6 


Both 


217 


93.3 


-0.70 


0.84 


19.6 


24.3 


4.62 


31.5 


No 


239 


93.3 


-0.70 


0.87 


37.6 


6.14 


4.10 


65.2 


No 


241 


91.5 


-0.68 


0.83 


32.1 


17.8 


7.20 


76.9 


No 


242 


91.5 


-0.69 


0.81 


19.7 


20.6 


15.0 


71.7 


No 


249 


87.1 


-0.63 


0.76 


10.0 


29.6 


4.87 


32.6 


Both 


257 


86.4 


-0.64 


0.75 


16.1 


15.2 


4.44 


61.7 


Yes 


258 


90.8 


-0.68 


0.82 


30.4 


21.9 


4.63 


54.1 


No 


260 


90.6 


-0.67 


0.82 


33.3 


14.1 


5.27 


42.8 


Yes 


267 


91.7 


-0.70 


0.83 


42.6 


19.6 


5.02 


68.4 


No 


268 


89.4 


-0.67 


0.80 


23.1 


12.7 


5.06 


24.4 


Yes 



Table 1: List of some selected orbits. Orbit number and wave properties. Prom left to 
right; degree of polarization, ellipticity, coherence, A2/A3, 9kB, |B|, Cone Angle av,B, 
Connection to the Martian Bow Shock (Yes, No, Both cases). 
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